Phase diagram of a Bose gas near a wide Feshbach resonance 
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In this paper, we study the phase diagram of a homogeneous Bose gas with a repulsive interaction 
near a wide Feshbach resonance at zero temperature. The Bose-Einstein-condensation (BEC) state 
of atoms is a metastable state. When the scattering length a exceeds a critical value depending on 
the atom density n, na 3 > 0.035, the molecular excitation energy is imaginary and the atomic BEC 
state is dynamically unstable against molecule formation. The BEC state of diatomic molecules 
has lower energy, where the atomic excitation is gapped and the molecular excitation is gapless. 
However when the scattering length is above another critical value, na 3 > 0.0164, the molecular 
BEC state becomes a unstable coherent mixture of atoms and molecules. In both BEC states, the 
binding energy of diatomic molecules is reduced due to the many-body effect. 
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I. INTRODUCTION 



Bose gases with Feshbach resonances have shown very interesting properties. The par- 
ticle loss rate due to the three-body recombination increases enormously near a Feshbach 
resonance A sudden change in the magnetic field generates oscillations between atoms 
and diatomic molecules Diatomic molecules can be produced when the magnetic field is 
either tuned through the resonance [3j or oscillates with a frequency corresponding to the 
molecular binding energy J4J]. Recently observed Effimov effect 5] may help to control the 
particle loss for more extensive studies near a Feshbach resonance. 

The atomic BEC state with strong interactions was explored in a variational approach 
6|. At resonance the energy per atom was found to be proportional to n 2 / 3 6j. A transition 
between atomic and molecular BEC states at the resonance was proposed 7|, [8|. However, 
these states were found unstable in some regions due to negative compressibility js]. Re- 
cently, a transition between the molecular BEC state and a coherent mixture state of atoms 

□ 

and molecules was proposed 10]. 

In a Bose gas with a Feshbach resonance, scattering states in the open channel are coupled 
to bound states in the closed channel. Eigenstates of diatomic molecules are superpositions 
of both open and closed channel states. Near a wide resonance where the effective range 
of the interaction is very small, the magnitude of the close-channel component is much less 
than that of the open-channel component, which can be seen in the renormalization of the 
molecular propagator [llj]. Thus an effective model of open-channel atoms should be able 
to describe the system with a wide resonance. A single- channel model has been widely used 



to describe both normal and BEC states of Bose gases 



13], given by 



fi 2 

H = ^V 2 ip + -^^W, (1) 

2/77; 2 

where ip is atomic field operator, g is the coupling constant, g = 4irh 2 a/m, and a is the 
scattering length. When the scattering length is negative, a < 0, the system is unstable and 
subject to mechanical collapse at low temperatures 14, Q|- Therefore in the following 



we focus on the region with repulsive interactions, a > 0. 

A crucial question about the single-channel model is whether or not it is capable of 
describing diatomic molecules near the resonance. The answer of this question is positive, 
because the Hamiltonian given by Eq. (jTJ) has bound eigenstates of two atoms at a > 0. 
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FIG. 1: The phase diagram of a Bose gas near a wide Feshbach resonance at zero temperature. 
The atomic BEC state (ABEC) is metastable with < \i < 0.52E a . The molecular BEC (MBEC) 
state is with -0.5E a < /x < —0.2lE a . 



These bound states are the eigenstates of diatomic molecules, given by 



1 



Af 



>|o), 



(2) 



VFV^ + 2e/ q+p/2 ^ q+p/ 

where V is the volume, AT is a normalization constant, e q = h 2 q 2 / (2m) , and the molecular 
binding energy in vacuum is given by E a = h 2 /(ma 2 ). The energy eigenvalue of the molecule 
is given by — E a + e p / 2, consistent with the energy of shallow bound states in the quantum 



scattering theory 



17|. It is important to note that the single-channel model is an effective 



field-theory model and ultravio 
correct renormalization scheme 



et divergences in calculations need to be removed by the 



3. 



In this paper, the phase diagram of a homogeneous Bose gas near a wide Feshbach 
resonance is obtained at zero temperature, as shown in Fig. [TJ The atomic BEC state is 
a metastable state when the chemical potential /i is between and 0.52E a , corresponding 
to < na 3 < 0.035. The atomic BEC state is unstable against molecule formation when 
na 3 > 0.035. The molecular BEC state is stable when —0.21E a > fi > -0.5E a , or 0.0164 > 
na 3 > 0. There is no solution when — 0.21E a < fi < 0. The coherent mixture state of atoms 
and molecules is subject to mechanical collapse, although // can be smaller than —0.5E a at 
very high densities. In both BEC states, the molecular binding energy is smaller than E a 
due to interaction. This phase diagram is valid when a is much larger than the effective 
range. 
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II. THE ATOMIC BEC STATE 



Bose atoms condense below the BEC transition temperature. The atomic BEC state in 
the dilute case is well described by Bogoliubov's theory in which the atomic field operator 
has a finite expectation value, (ip) = ipQ. In the grand canonical ensemble, the grand 
thermodynamic potential F = H — yL^tj) has a constant part Fq = gn^/2 — firio, where /i is 
the chemical potential and n = \ipo\ 2 . The saddle-point condition 5F /5ijj yields \i = gn . 
For convenience ipo can be chosen to be positive. The quadratic part of the grand potential 
describes the gaussian fluctuation in the atomic field 5tp = ip — tpo, 

h 2 1 

F 2 = 5^V 2 5i) + gn 5^Sip + -gn (5^5^ + h.c). (3) 

2m 2 

By Bogoliubov transformation, the quadratic grand potential F 2 can be diagonalized, 

F 2 = C + 5> fc 4c k , (4) 

k 

where the quasi-particle energy is given by E k = ^Je^k + 2gn ), = h 2 k 2 /(2m), the 
quasi-particle operator is given by Ck = Ukipk — Vkip^, u\ = [1 + (e^ + gno)/Ek}/2, and 
v k = u l ~ 1- The quantum depletion, the vacuum of quasi-particles, contributes to the 
ground-state energy by C = 8gno(mgno/h 2 ) 3 ^ 2 / '( 157r 2 ) . The higher order terms in the grand 
potential, 

g8^8ip[(ij) 5ij) + h.c) + 5^5^/21 

n 

are neglected in Bogoliubov's approximation. In Popov's approximation [12J , the mean- field 
terms of non-condensed atoms are added in the particle-hole channel, and the chemical 
potential is shifted, fi = g(n + 25n), where Sn = 8n ^Jn a 3 /ir/3 is the the atom density of 
the quantum depletion. In the diagrammatic approach in the traditional theory of a dilute 
Bose gas [ljj, thermodynamical properties can be calculated by perturbation in the order 
of beyond Bogoliubov's approximation. 

The traditional theory is accurate in the dilute case, where the perturbation can be 
stopped at a sufficient order of yna? without losing much accuracy. However near the reso- 
nance, there are strong fluctuations around the condensate. In addition to single-atom exci- 
tations, collective excitations must be considered, such as excitations of diatomic molecules. 
The dispersion of molecular excitations can be obtained from the poles of the two-particle 
correlation function of non-condensed atoms. In the atomic BEC state, the particle-particle 
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channel and the particle-hole channel are coupled, and the correlation function is a 3 x 3 
matrix given by 

XaP (v -v',t- t') = ~{T[b a (r,t)bl(r',t')}), (5) 

where b\ = ip 2 , 62 = b\, 63 = 2ip^ip, and T is the time-ordering operator. 

Due to the strong interaction near the resonance, any perturbation theory truncated at 
any finite order of \/na 3 will fail. In the following, we adopt the Random Phase Approxi- 
mation (RPA) which focuses on the renormalization to the two-body interaction. In RPA, 
there are Feynman diagrams from all orders of perturbation and the result can often be 
applied to the strong- interaction region with correct qualitative features [131 ]. 

In RPA, the correlation function is given by 

X (k,^) = [l-^ )(k,^)]-V )(k,^). (6) 

The function %(°)(k, to) is the correlation function calculated in Bogoliubov's approximation, 
given by 

v (o) rk n] _ f d 3 k' ul,uf k _ kll vlv%^\ 
Xn {KUJ) ~ J (2^ ( A(k,k',c) + A(k,k',-co) h 

Xl2{K") - 2 J (2?r) 3 B (k,k',u) ' 
(0), _ f d 3 k' Mfc'ffe'Mji_ k '| Uk'Vk'vfk-k'l 



(2vr) 3V A{k,k',u) A{k,k',-u) 



v (0), t \ f d3k ' K^|k-k'| +^'M|k-k'|) 2 (7 , 

X33 {k,uj) - 2 J — B (k,k',u) ' [) 

X&fcw) = X11 (-k, -u), X2?( k >^) = x2 (- k > and xi° } (k,w) = xfifav) for 3 ^ *» 
where A(k,k',uj) = Tiu - E k - E\^ w \+i8 and 1/B(k,k',u) = 1/A(k, k', u) + 1/A(k, k', -co). 

The dispersion of diatomic molecules can be obtained from the pole of the correlation 
function given in Eq. 

det|I-£ X (0) (k,^)| = 0, (8) 

where I is the identity matrix. The real part of the molecular excitation energy at k = 
is the negative of the binding energy which is plotted in Fig. [2j When the interaction 
is weak, na 3 <C 1, the molecular binding energy is approximately given by E a — 2gno, 
which recovers the vacuum result in the dilute limit. The size of the bound state is of the 
order of the scattering length a. When the wavevector k is of the order of 1/a or larger, 
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FIG. 2: The binding energy and damping of diatomic molecules versus n 1 / 3 a in atomic and molec- 
ular BEC states. The solid line is the binding energy and the dash-dotted line is the imaginary or 
damping part of the molecular excitation in the atomic BEC state. The dashed line is the binding 
energy in the molecular BEC state. 

the atomic-excitation energy is approximately given by e\. + gn . This energy shift gn in 
atomic excitation energy at high momentum causes the reduction in the molecular binding 
energy by -2gn . 

When the interaction is stronger, the imaginary part of the molecular excitation energy 
becomes noneligible. At the critical value of na 3 = 0.035 or // = 0.52E a , the real part of 
the molecular energy vanishes and only the imaginary part is finite, implying that beyond 
this point the system is totally unstable against molecule formation. This instability is a 
many-body effect, due to the process described in RPA in which two atoms from two pairs in 
the quantum depletion form a molecule leaving the rest two atoms in the two pairs excited. 
Such process is absent in the normal state where there is no off-diagonal correlation. In 
contrast, the particle loss due to the three-body recombination is a few-body effect which is 
present in both BEC and normal states. 

III. THE MOLECULAR BEC STATE 

When the scattering length is positive, a > 0, diatomic molecules have lower energy than 
atoms, and can condense into a molecular BEC state at low temperatures. In the single- 
channel model, the molecular condensation can be described by the off-diagonal long range 
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order A = g(ipip) 7^ 0. For simplicity, we assume A > 0. The mean- field grand potential 
describing the molecular BEC state is given by 

F m = - |VvV + (2gn - + y (V>¥ + W> ( 9 ) 

where i 7 ^ = — {gn 2 + A 2 /2g). This grand potential _F m can be diagonalized by Bogoliubov 
transformation, 

F m = F mQ + E Efccj.Ck, (10) 
k 



where E k = \J{e k + 2gn - /i) 2 - A 2 , c k = u k ipk - v k i(>l_ k , u\ = [1 + (e fc + 2#n - n)/E k }/2, 
v k — u 2 . — 1, and the constant F m0 is the ground state energy 



1 A 2 



k 

The parameter A can be determined self-consistently, 

d 3 k 



A = g J u k v k 



or 

1 r d 3 k r 1 1 



y 7 (2vr)3 L 2£ fc + 2eJ' (U) 
where the last term in the integrand l/(2e k ) is a counter term. The chemical potential \i 
and parameter A can be solved from Eq. (jlip and the equation for the density n, 



n 



d 3 k 2 /■ <i 3 £; e fc + 2gn - [i 

(2^ Vk = J j2^ [ —2E k 1] " (12) 



There is a gap in the atomic excitation given by Eq = J {2gn — /1) 2 — A 2 , consistent with 
the fact that atoms have higher energy than molecules. The molecular binding energy in 
the molecular BEC state is given by 2Eq, equal to the energy difference between two atomic 
excitations and a molecular excitation at k = 0. The binding energy is plotted in Fig. [21 In 
the dilute limit where na 3 <C 1, the gap E is approximately equal to E a /2. However at the 
critical value of na 3 = 0.0164, the gap vanishes. When na 3 > 0.0164, there is no solution, 
indicating that the molecular BEC state no longer exists. 

Similar to the case in the atomic BEC state, the dispersion of molecular excitations in the 
molecular BEC state can be obtained from the poles of the two-particle correlation function 
X- In RPA, the pole is given by Eq. ([8]), and the mean- field correlation function x is given 
by Eq. ©, with the coefficients, u k and v k , and the excitation energy E k replaced with 
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FIG. 3: The molecular velocity v m versus n 1//3 a in the molecular BEC state. At na 3 = 0.0164, the 
velocity v m vanishes. 



the values in the molecular BEC state. When k = and uj = 0, Eq. (jSj) is automatically 
satisfied following Eq. (J7J [TT]) . showing that the molecular excitation is gapless. At small 
k and u, to the leading order of k, the molecular excitation frequency is linearly dispersed, 
ojk ~ v m k, where v m is the molecule velocity. 

The molecule velocity is plotted in Fig. [3j In the dilute limit when na 3 <C 1, the molecule 
velocity is approximately given by 



h 



m 



(13) 



If Eq. ([TBI is compared to the phonon velocity v p = h^Anna/m in the atomic BEC state, 
the naive estimation of the molecule-molecule scattering length is 6a, larger than the result 
from solving the four-body problem [lj]]. This reason for this discrepancy might be that in 
the dilute limit the mean-field density given by Eq. (fl~2l) is smaller than the true density due 
to fluctuations, similar to that in the molecular BEC state of the Fermi gas |20| . When the 
interaction is stronger, at the critical value of na 3 = 0.0164 or fi = —0.21E a , the molecule 
velocity v m drops to zero, which is very similar to softening of phonon modes in solids near 
structural phase transitions. The vanish of both the atomic excitation gap Eq and molecule 
velocity v m at na 3 = 0.0164 indicates that beyond this critical point the molecular BEC 
state does no exist. 
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IV. THE COHERENT MIXTURE OF ATOMS AND MOLECULES 



In both BEC states, the molecular binding energy is reduced by the many-body inter- 
action. In the region beyond these BEC states, the system may contain both atoms and 
molecules. Here we consider a coherent mixture state of the atomic condensation with 
■ip = (ip), and molecular condensation with A = g(5ip5ij)}, where Sip = ip — ipo. 111 this 
mixture state, the two parameters ipo and A are in principle independent, whereas in the 
atomic BEC state the parameter A is a function of ipo. 

In the mean-field approximation, the grand potential is given by 

F mx = FW - ^S^V 2 SiP + {2gn - fidtffy 

+I[(A + # 2 )#W + M, (14) 

where Fj® = — [g5n 2 +n {[i— gn /2) + \A\ 2 / (2g)], n = \ip \ 2 , Sn = (Sip^Sip), and n = n +5n. 
Its mean-field expectation value, F mx = (F mx ), is given by 

F m x = | |-0o + ^"l 2 + 2gn 5n + |tfn 2 - fi(5n + n ). (15) 

The parameter ipo should minimize the grand potential, dF mx /dipl = 0, which yields the 
saddle-point equation 

fx = g(n + 25n) + A^. (16) 

Wo 

Apparently, the product Aip^ 2 is real. For simplicity, we choose A to be positive, which also 
means that i\) 2 is real. 

The mean-field grand potential can be diagonalized by Bogoliubov transformation, where 
the field operator of the quasi-particles are given by Ck = Ukipk — v k ip_ k , the coefficients are 
given by u\ = [1 + (e fe + 2gn — fi)/E k ]/2 and v\ = u\ — 1, and the quasi-particle energy is 
given by E k = ^ {e k + 2gn — /i) 2 — (A + gipo) 2 - The parameter A can be determined from 
the self-consistency equation, 

d 3 k 



i. e. 



f « fc 
A = 9 J j2nf UkVk > 

ifj 2 f d 3 k 



g A + # 2 J (2tt) 
The expression of the total density is now given by 



( r 1 1 i (m\ 



r d 3 k 1 e k + 2gn-n 
n = no + i (2^2 ( E k 1} - (18) 
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FIG. 4: The chemical potential // versus n 1//3 a. The solid line is ji in the atomic BEC state in 



121 ] for na 3 < 0.035. The dashed line is /i in the molecular BEC state 



Popov's approximation 
which exists at na 3 < 0.0164. The dotted line is [i in the mixture state. 

For a fixed total density n, both A and ipo can be solved from Eq. (fPo| [TTI [TS]) . In the 
solution, A and have opposite signs, whereas in the atomic BEC state A and ip^ have 
the same sign. At the transition point to the molecular BEC state, the atomic condensation 
density no is zero. The solution shows a reentrant behavior in a tiny region with 0.0156 < 
na 3 < 0.0164. When na 3 > 0.0164, the solution is unique. However, in this region as shown 
in Fig. HI the chemical potential /i decreases with the increase in density n, which means that 
the compressibility dfi/dn is negative and the mixture state is subject to mechanical collapse. 
This instability is very similar to the instability of Bose gases with attractive interactions 
, Q|- In a trap, the mixture state may be stabilized by the finite-size effect under 
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certain conditions, as discovered in trapped Bose gases with attractive interactions 21]. 

Conclusion and Discussions. The phase diagram of a homogeneous Bose gas near 
a wide Feshbach resonance is studied at zero temperature. In the atomic BEC state, the 
real part of the molecular excitation energy vanish at na 3 = 0.035, implying the instability 
against molecular formation for stronger interactions. In the molecular BEC state, the 
atomic excitation energy is gapped and the molecular excitation energy is linearly dispersed. 
Both the gap and molecular velocity vanish at na 3 = 0.0164, above which the molecular 
BEC state no longer exists. In both BEC states molecular binding energies are reduced 
by the many-body interaction. In the coherent mixture state of atoms and molecules, at 
na 3 > 0.0164, the compressibility is found to be negative, indicating that the mixture is 
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subject to mechanical collapse, similar to Bose gases with attractive interactions. 

It is an open question whether or not there are more exotic states in regions beyond the 
BEC states. So far three-body interactions are ignored in the BEC states. Although they are 
crucial to the dynamical properties such as the particle-loss rate and the stability time, the 
quasi-equilibrium properties is unlikely affected. However triatomic Effimov molecules may 
form near the resonance due to three-body interactions. Whether or not a gas of Effimov 
triatomic molecules can exist is at present unknown. 
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